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We study the 0(n) loop model on a dynamically triangulated disk, with a new type of 
boundary conditions, discovered recently by Jacobsen and Saleur. The partition function 
^h. of the model is that of a gas of self and mutually avoiding loops covering the disk. The 
Jacobsen-Saleur (JS) boundary condition prescribes that the loops that do not touch the 
boundary have fugacity n G [—2, 2], while the loops touching at least once the boundary are 
given different fugacity y. The class of JS boundary conditions, labeled by the real number 
y, contains the Neumann (y = n) and Dirichlet (y = 1) boundary conditions as particular 
cases. Here we consider the dense phase of the loop gas, where we compute the two-point 
boundary correlators of the L-\eg operators with mixed Neumann-JS boundary condition. 
The result coincides with the boundary two-point function in Liouville theory, derived by 
Fateev, Zamolodchikov and Zamolodchikov. The Liouville charge of the boundary opera- 
tors match, by the KPZ correspondence, with the L-leg boundary exponents conjectured 
by JS. 
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1. Introduction 



The solvable statistical models that have a geometrical description in terms of self and 
mutually avoiding clusters, like Ising, 0(n) and Potts models, can be also formulated 
and solved on a dynamical lattice @ @ [§J || • A statistical model defined on a 
dynamical lattice is said to be coupled to gravity, since the sum over lattices gives a 
discretization of the path integral over Riemann metrics on the world sheet. For each 
critical point described by a 'matter' CFT, the 'coupling to gravity' consists in adding a 
Liouville and ghost sectors and dressing the scaling operators by exponents of the Liouville 
field 0@ |Ty]. The description of the critical points based on Liouville theory allows to 



interpret the wealth of exact results about statistical models on dynamical lattices obtained 
via matrix-model or combinatorial techniques. 

A statistical system on a random lattice exhibits qualitatively the same critical phases 
as that on a regular lattice, but with different critical exponents. At a critical point 
characterized by a conformal anomaly c < 1, the conformal weight h of a matter fieldS can 
be extracted from its 'gravitational dimension' A, which determines the scaling properties 
of the correlation functions involving this field |L(J , 



A A - 7str ) 7s 2 tr 
h = — , c = 1 - 6- . (1.1) 
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The exponent 7 str (gamma-string) describes the critical fluctuations of the area of the 
random surface. 

The solutions of the KPZ scaling relation ( |1 . 1| ) can be parametrized by a pair of 
numbers r and s, not necessarily integers: 

_ ( r/6 _ g6) 2_ (1/6 _ 6) 2 _ r/b _ sb _ {1/b _ b) 

"rs i^rs , 5 t-^rs nl ■ \ J 

4 2b 

The dependence on the matter central charge is through the positive parameter b < 1, 
defined as 

b 2 = — ^ < 1. (1.3) 

1 - 7str 

The conformal weights have the symmetry h rs = h- r> - s , unlike the gravitational dimen- 
sions. To each conformal weight one can associate two gravitational dimensions, A rs and 
A_ rj _ s , which are the two roots of the quadratic relation (|1.1| ). They correspond to the two 
possible gravitational dressings of the matter conformal field by Liouville vertex operators. 

The correspondence between the critical phenomena on flat and dynamical lattices is 
particularly useful in presence of boundaries. The boundary critical exponents O] on flat 



Only operators without spin (h = h) survive after coupling to gravity. 
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and dynamical lattices are again related by ( Jl.lp . Boundary 2D quantum gravity became a 
powerful tool for evaluating exact critical exponents complementary to the Coulomb 
gas techniques [[T3| and the Schramm-Loewner Evolution (SLE) For systems coupled 
to gravity it is possible, by cutting open the path integral along cluster boundaries, to 
solve analytically problems which whose exact solution is inaccessible on flat lattice. On 
the other hand, our understanding of the boundary phenomena in 2D gravity is much 
helped by the progress achieved in boundary Liouville theory in the last years [IS] |TB[ [TTl . 



In this paper we will focus on the possible boundary conditions and the spectrum of 
boundary exponents of the O(n) model coupled to 2D gravity 0. Our main purpose is to 
check a very interesting conjecture made in |18[ about the general boundary conditions for 
the 0(n) model on a flat lattice. 

The O(n) model has a continuum transition if the number of flavors is in the interval 
— 2 < n < 2, where it can be parametrized by an angle, 

n = 2cos(7r#). (1.4) 

The boundary 0{n) model was originally considered with Neumann boundary condition 
(the loops are reflected from the boundary). The boundary scaling dimensions of the L-leg 
operators, realized as sources of L open lines, were conjectured in [TJ] and then derived in 



. Furthermore, the partition function of the 0{n) model on the annulus with Neumann 



boundary conditions was evaluated in |21[|. Another obvious boundary condition is the 



Dirichlet boundary condition, studied in [22] []23j , for which there is an open line ending at 
each site of the boundary.i The dimensions of the L-leg boundary operators with mixed 



Dirichlet and Neumann boundary conditions were computed in ||24j| |[23|| by coupling the 
model to 2D gravity and then using the KPZ scaling relation ( |1 . 1| ) . 

Recently, Jacobsen and Saleur put forward a proposal about the complete classification 
of the boundary conditions of the 0(n) loop gas model in the dense phase, described by, 
in general non-rational, CFT with b 2 = 1 — 9. Their proposal is based on a previous 
work [ 25 1 and possibly overlaps, for rational 9, with the results of |2jJ. The claim of [|18 



is that there is a continuum of boundary conditions characterized by a real variable y. 
The Jacobsen-Saleur boundary condition, which we will denote in the following by JS, is 
defined by counting the loops that touch at least once the boundary with different fugacity 
y, while the loops that do not touch the boundary are counted with fugacity n. The JS 
boundary conditions contain as particular cases Neumann [y = n) and Dirichlet (y = 1) 
boundary conditions for the 0(n) field. 



2 In these papers the loop gas was considered in the context of the SOS model, for which the 
Dirichlet and Neumann boundary conditions have the opposite meaning. 
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In order to evaluate the L-\eg exponents, the authors of |TJ[] considered the loop gas on 
an annulus, with Neumann boundary condition on the inner rim, JS boundary condition 
on the outer rim, and L non-contractible loops separating the two boundaries. The loop 
model with these boundary conditions were called there boundary loop model (BLM). 
According to ||18|| , BLM with L non-contractible lines has two sectors, the blobbed and 
the unblobbed one, characterized by two different for (L > 1) critical exponents. In the 
blobbed (unblobbed) sector the outmost non-contractible line touches at least once (does 



not touch) the outer rim. It is argued in [T^] that the scaling exponents in the two sectors 



are characterized by the conformal weights 

unblob 

h L = /i_ r) _ r _|_£ 



blob 



(dense phase), (1-5) 



where the real parameter r is related to the boundary fugacity y byi 



, , sinlYr + 1)tv9] 
sm(mv) 



The boundary exponents for Neumann-Neumann |H| and Neumann-Dirichlet bound- 
ary conditions appear as particular cases of ( |1.5| ) when r = 1 and r = correspondingly. 

Inspired by |f[8 |, in this paper we analyze the JS boundary conditions for the 0{n) 



loop gas coupled to 2D gravity. Put in string theory terms, the problem we address is to 
describe the ensemble of D-branes in bosonic string theory whose target space is the (n — 1)- 
dimensional sphere. We restrict ourselves to the dense phase of the 0(n) model, where we 
derive loop equations in the form of recurrence relations between the the boundary two- 
point functions of the L-\eg and (L — l)-leg operators. The loop equations are obtained 
by cutting open the world sheet along segments of loops that connect two points of the 
boundary. In the continuum limit, the solution of the loop equations is given, up to a 
normalization, by the two-point boundary correlator in boundary Liouville theory fl5 



The scaling exponents characterizing the solution reproduce, via KPZ relation (|1 . 1|) , the 



spectrum of boundary conformal weights conjectured by JS [jl8 |. 

The paper is organized as follows. In sect. 2 we give a short review of the results we 
will need about the loop gas coupled to gravity. In particular, we derive the loop equations 
for the disk partition function with Neumann boundary conditions. In sec. 3 we formulate 
the JS boundary conditions in terms of the boundary measure for the 0(n) spins and 
derive the loop equations for the two-point correlators with mixed Neumann-JS boundary 
conditions. In section 4 we take the continuum limit of the loop equations and evaluate 
the boundary L-leg exponents. Here we also show that the loop equations can be cast in 
the form of the functional equations, generated by the boundary ground ring of Liouville 
gravity. Summary of the results and some concluding remarks are presented in Section 5. 



Our notations are related to those used [IS] by = 2 — 1 



p+i 
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2. The boundary 0(n) loop model coupled to 2D gravity 



2.1. Disk partition function with Neumann boundary conditions 



The 0(n) model |L3|, defined originally on a regular hexagonal lattice, can be considered 
on an arbitrary trivalent planar graph T*, which is dual to a triangulation of the disk T. 
The local fluctuating field is a n-component spin S(r) associated with the vertices r of T*. 
The Boltzmann weight for given configuration is 

J] (t + S(r) ■ S(r')) , (2.1) 

<rr'> 

where the product is over the links <rr'> of T*. Here we consider the maximally packed 
dense phase t = 0. The partition function for given triangulation is the trace over these 
weights, defined by Tr[l] = 1, Tr[S a (r)S b (r)] = 5 ab and Tr[S a {r)] = Tr[S a (r)S b (r)S c (r)} = 
0. Expanding the trace as a sum of monomials, the trace for given triangulation T can be 
written as a sum over all configurations of densely packed self-avoiding, mutually avoiding 
loops on T* (Fig.l). Each loop is taken with a weight n. 



Fig. 1: Loops on a trivalent planar graph T* (left) dual to a triangulation T (right) 



The disk partition function of the 0(n) model coupled to 2D gravity is defined as the 
double average: with respect to the 0(n) field on given triangulation and with respect 
a sufficiently large class of triangulations. We will consider all possible triangulations, 
including the degenerate ones, which are dual to trivalent planar graphs with several 
connected components. In such a triangulation a boundary edge can either belong to a 
triangle, or be identified with another edge of the boundary. The measure in the ensemble 
of triangulations of the disk is determined by the 'cosmo logical constant' x coupled to 
the area and the 'boundary cosmo logical constant' ( coupled to the boundary length. By 
definition the area of a triangulation is the number of its triangles and the boundary length 
is the number of edges along the boundary. 
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The basic observable in the 0(n) model coupled to 2D gravity is the disk partition 
function x). Its derivative W = —d^ is given by the series 



W(C) = J2C l Wi(x), (2.2) 

where Wi is the non-normalized expectation value of the Boltzmann factor (|2.1|) in the 
ensemble of triangulations % with boundary of length I. Since there is no restriction 
for the 0(n) spins at the boundary, we have Neumann boundary condition. The loop 
expansion of W\ is 

— x ~(Area of 71) ^(Number of loops) ^ 3^ 

71 loops on T* 

Note that the here boundary has a marked point, hence there is no symmetry factor l/l m 
the sum. The generating function ( |2.2[ ) is the boundary one-point function of the identity 
operator. 



2.2. Loop equation for the disk partition function 

Here we remind the combinatorial derivation of the loop equations for the the disk am- 
plitude with Neumann boundary conditions ( |2.2| ). This derivation, first given in ||, is 
a useful exercise to do before passing to the more complicated case of mixed boundary 
conditions. 

The triangulations filled by loops that enter in the sum in the r.h.s. of (|2.3| ) can be 
divided into two classes (Fig. 2). The first class comprises the degenerate triangulations 
for which first boundary edge is not an edge of a triangle, but is connected directly to the 
/c + l-th boundary edge. The contribution of such triangulations factorizes to Wk-iWi-k-i- 
There are I — 1 such terms, with k = 1, 2, / — 1. For the rest of the triangulations entering 
the sum in ( |2.3| ) the first edge belong to a triangle, which must be is visited by a loop, 
since the loops are densely packed. Now consider the ensemble of all triangles visited by 
this loop. These triangles form a closed strip that contains the loop. Let q and p be 
respectively the lengths of the internal and external boundaries of the strip. For given q 
and p the expectation value factorized to the contribution of the internal disk W q , that of 
the external disk Wi-i- p , and the number of realizations of the strip, y . We have also 
a factor n because of the loop contained in the strip. Adding the contributions of the two 
classes of triangulations, we obtain the following bilinear equation for the Wk, 

k=l P,q>0 P'^' 
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where by definition Wo = 1. Equations of this type are usually called loop equations, in 
terms of the generating function ( |2.2| ), the loop equation reads 

[-CW(0 + W(C) 2 + nW(x - C)W(0] < = 0, (2.5) 
where [ ] < denotes the negative part of the Laurent series in (. 




Fig. 2: Loop equation for the disk partition function with Neumann boundary 
conditions. The small rectangle stands for the marked point on the boundary. 



It is consistent to assume that the solution has a cut [a, b] on the real axis, with a < b < x/2. 
Then one can write the projection [ ]< in (|2.5| ) as a contour integral, 

1 - iW (0 +W H0 + n<ff m0^m w{x _ = o. (2.6) 

J Ztvi z — z 

where the contour of integration encircles the cut of W(() and leaves outside the cut of 
W(x — C). Equation ( |2.6|) yields a condition for the discontinuity across the real axis 



DiscW(C) hC + W(C + i0) + W(C - i0) + n W(x - ()] = 0, (( G 1). (2.7) 

The condition ( |2.7p , after being symmetrized with respect to £ — > x — C, implies that 
certain bilinear combination of VL(C) and W(x — Q has zero discontinuity on the real axis 
and therefore is analytic in the whole complex plane. This leads, taking into account the 
asymptotics W(C) ~ 1/Q at infinity, to the functional identity 0] 

W(C) 2 + W(x - C) 2 + nW(()W(x - C) = (W(() + {x- QW(x - () - 2. (2.8) 



3. Boundary correlators with mixed Neumann/ JS boundary conditions 

3.1. The JS boundary conditions in terms of spins and loops 

The JS boundary condition can be introduced by restricting the 0{n) spins on the bound- 
ary to take values in a submanifold of 'dimension' y. That is, the first y components the 
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0(n) are free, while the rest n — y components are fixed. This is equivalent to replacing 
the Boltzmann weight ( |2.1| ) with 



II S(r)-S(r') f[ J2S b (r)S b (r). 



(3.1) 



<rr'> 6 bulk 



rEboundary b=l 



The Boltzmann weight ( |3.1| ) is invariant with respect to a subgroup 0(y) C 0(n), which 
means that the JS boundary conditions are associated with the conjugacy classes of 0(n). 
This is a necessary condition to have conformal invariant boundary theory p7|. 

In the original formulation of the 0(n) model both n and y are integers, § but the 
result of evaluating the trace can be analytically continued for non-integer values of y and 
n. The disk partition function is then formulated as a gas of fully packed loops on the 
world sheet, having two different fugacities. The loops that do not touch the boundary 
have fugacity n, while the loops that are reflected from the boundary one or several times 
have different fugacity y. 

There are two classes of local operators compatible with the boundary measure in 
( |3.1| ), which we denote by and S]k They are defined as 0(y) invariant polynomials of 
the spin components: 



/] Set! ■ ■ ■ S, 



O.L 1 



l<a 1 <...aL<y 



if 



Sb 1 ■ ■ ■ Sb L 



y+l<6i<...bz,<r 



(3.2) 





Fig. 3: Topology of the loop configurations for the boundary two-point functions. 



We will investigate the boundary correlation functions of the operators ( |3.2| ) with 
Neumann boundary condition on the left segment and JS boundary conditions on the 
right segment of the boundary. We denote these correlation functions by 

DliU) = ( c N [4]i[4i V. k , DUtf) = (SlsilLPili ) H . ■ (3-3) 

' disk \ / disk 



4 It is always possible to consider part of the components of the O(n) vector as anticommuting 
variables, so that the restriction < y < n can be avoided. 
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The symbol ( ) disk means a double sum over the 0(n) spins and over the triangulations 
of the disk, characterized by the boundary cosmological constants ( and ( associated with 
the two segments of the boundary. In terms of the loop gas, the two-point functions ( |3.3| ) 
are the partition function of the loop gas with L open lines connecting two marked points 
on the boundary, as shown in Fig. 3. In the loop expansion of D L , the configurations 
where the rightmost open line touches the JS boundary have the same weight as those in 
which it doesn't. In the case of Sj;, the lines that touch the JS boundary have zero weight. 
In the terminology of ||18|| , the correlators D L and are the two-point functions of the 
L-leg operator respectively in the blobbed and unblobbed sectors. 



3.2. Loop equations for the L-leg correlators with mixed N/JS boundary conditions 
The correlation functions ( jOl) are defined as the series expansions 

oo 

D L (C, C) = E Dfi C 1 - 1 C ! -\ (3.4) 
i, i=o 

where Dl can be either or Dj^. The coefficients D^~ are the the two-point functions 

with fixed boundary lengths, / and /. The functions Dl and Dl-i satisfy a simple re- 
currence relation [E23j . This relation in fact holds for Neumann boundary condition on 
the left segment and any boundary condition on the right segment of the boundary. It 
is obtained by taking into account all possible configurations of the strip containing the 
leftmost open line (Fig. 4). 




Each such configuration is characterized by the lengths p and q respectively of the left 
and the right boundaries of the strip. The length of the leftmost line is then p + q. 
The strip splits the triangulated disk into two pieces in such a way that the Boltzmann 
weights factorize, so that the sum over triangulations and loops can be done in each piece 
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separately. The left piece contributes the disk amplitude Wi +P , while the right piece yields 



the boundary two-point function D - 1 , again with Neumann boundary condition on the 
left segment. Taking into account that the strip in the middle can be realized in (p+q)\/p\q\ 
ways as well as the the factor x~ p ~ q associated with the area p + q of the strip we get 

^r=E^ fi p#*- M ^r 1 - (3-5) 

p,q=0 F ' H ' 



Written in terms of the generating functions Q3.4Q , this equation takes the form 

D L (CO = [W(C)D L - 1 (x-CX)} < (3-6) 

where we used the same notations as in ( |2.5|) . We again express the projection [ ]< as a 
contour integral: 

This gives for the discontinuity on the real axis 

Disc CJ D L (C, C) + Disc c W(C) ■ D L _ x {x - C, C) = 0. (3.8) 



Equation ( |3.5| ) was derived for L > 2, but we will extend the definition of the two- 
point functions to L = 0, so that the it holds also for L = 1. The subtle point here is how 
to weight the degenerate triangulations where the left and right boundaries touch at one 
or several points. The two-point function Dq is defined as a sum over all triangulations, 
including those where two boundaries touch, while in the loop expansion of the two-point 
function Dq the triangulations with touchings are excluded. Such partition functions be- 
have as though local operators were inserted at the marked points O] . We will think of Dq 



and Dq as the boundary two-point function of the operators Sq and Sq, correspondingly. 

The recurrence relations ( |3.5|) allow to determine all correlation functions ( |3.3|) , once 
the functions 

A = Dq and B = D 1 } 

are known. Below we derive two independent non-linear equations for A and B and find 
their unique solution in the continuum limit. These equations involve nontrivially only 
the cosmological constant Q of the Neumann boundary, while the cosmological constant 
C of the JS boundary enters as a parameter. Therefore from now on the dependence on 
Q will be implicit in our notations. The first equation is derived by splitting the sum of 
triangulations contributing to A l j into four sets associated with the way the first edge of 
the JS-boundary is connected to the rest of the triangulation (Fig. 5). First, there is the 
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possibility that the JS boundary has zero length (Z = 0). In this case A { ,~ = Wi. If the 
JS boundary contains at least one edge, we count 2 possibilities. The first edge can be 
glued to another edge of the JS boundary, or it can be the edge of a triangle containing 
a segment of a loop. The last possibility is realized by two types of configurations: (a) a 
loop that touches the boundary only once and (b) a loop that touches the boundary at 
least twice. 




In the case (a) we get a product A l+p jW q , as in the loop equation (p.4|) for the disk 
amplitude. In the case (b) we can apply the same argument as the one used in the 
derivation of ( |3.5| ). The strip that splits the disk into two here consists of the triangles 
visited by the most external arc of the loop. The left piece is bounded by the Neumann 
boundary, the left boundary of the strip of length p and a piece of the JS boundary of 
length k. It yields a factor Ai +p ^- The right piece is the partition function of the disk with 
mixed Neumann/ JS boundary conditions and an open line connecting the extremities of 
the two segments of the boundary. The open line can touch both boundaries unrestricted 
number of times. By definition this is the two-point function B q j_ k _ 2 , where q is the 
length of the right side of the strip. 
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Summing up the four terms and taking into account the combinatorial factors, we get 

1-2 



k=0 



+ E ^-^WA^r-i^ (3.9) 

p,q = 

1-2 oo 



-k-2- 

k=0p,q=0 



The equation satisfied by B is very similar to ( |3.5| ) . Here we distinguish two possibil- 
ities: the open line does not touch (touches at least once) the JS boundary (Fig. 6): 



p,q=0 F H 



(3.10) 

y> (p + g)! x -p- q -i B . A h 

+ X n i+p,l-k-l A q,k- 



plql 




In terms of the generating functions ( |3.4| ) the loop equations ( |3T9| ) and ( |3.10|) state 

C~A(C) = w(0 + A(c)w(C) + y [A(0 (W(x -() + b(x - c))] < (3.n) 

fl(C) = [^(C)A(x - C)]< + [B(C)A(x - C)]< • (3.12) 

Representing, as before, the projection [ ]< as a contour integral and taking the disconti- 
nuity across the real axis, we obtain 

2/Disc c A(C) • (W(x - C) + B{x - C) + ±(W(C) ~ C)) + Disc c W(C) = (3.13) 

Disc c S(C) = Disc c W(C) • A(x - C) + Disc c S(C) • A(x - ()• (3.14) 
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3. 3. Final form of the loop equations for D'j and Dq 

The linear terms in equations (|3.13|) and (|3.14|) can be eliminated by shifting the observ- 
ables A = Dq and B = D±. We redefine Dq and as 

D l }(():=B(() + W(0 + -(W(0-0, 

y (3.15) 

So (0 :=^(C)-1- 

The terms added only change the weight when one or both boundaries degenerate to a 
point and thus does not affect the critical behavior. For the shifted quantities ( ^.15| ), 
equations (|3.13j ) and ( |3.14| ) simplify to 

Disced (C) • D 1 } (x - C) + ^Disc C W(Q = , (3. 16) 

Disc c D^(z) ■D L (x-() + Disc c W(C) = 0. (3.17) 

The set of equations (|3.8|), ( |3.1fci| ) and ( |3.17| ) is overdetermined. Equation ( |3.17j ) is 
consistent with ( |3.8| ) for L = 1 under the condition 

Dj = ~. (3.18) 

This condition has simple geometrical explanation. Indeed, it is easy to see that the 
correlation functions Dq and Dq are, unlike their counterparts on a flat lattice, different. 
The difference comes from the fact that on a dynamical world sheet the Neumann boundary 
can come close to the JS boundary and touch it one or several times. Microscopically, the 
algebraic relation ( |3.18| ) can be understood as a geometric progression obtained by taking 
into account all possible touchings, 

D l Q l = 1 + A + A 2 + ... = — !— . (3.19) 



The two equations ( |3.16[ ) and ( ^.17| ), together with the asymptotics at infinity that 



follows from the expansion (|3.4| ), imply the following functional identity: 

D$(C)Dfr{x - C) = -W(C) ~ ~W(x - C) - -(W(0 - 0- (3-20) 

y y 

The loop equations ( |3.8|) , ( |3.16| ) and ( p.!7|) obtained here allow in principle to compute 
all boundary two-point correlators with mixed Neumann/JS boundary conditions. The 
three couplings, x, ( and (, associated respectively with the area of the triangulation, the 
length of the Neumann boundary and the length of the JS boundary, enter in the loop 
equations implicitly through the disk amplitude W(Q. 
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4. The continuum limit 



In this section we will study the continuum limit of the solution, in which the three cou- 
plings are tuned close to their critical values, £*, and £*. The solution in the continuum 
limit depends on the three renormalized couplings 

/Lt ~ x — x*, — (* and — C*- 

^.i. .DzsA; one-point function with Neumann boundary conditions 
The power series (|2.3| ) converges for x > x* , where 



x* = 2^2(2 + n) 

is the the critical value of the cosmological constant J7|. We introduce a small cutoff 
parameter (elementary length) a and define the renormalized cosmological constant fj,, 
boundary cosmological constant z and loop amplitude w as followsi 

_o 8 x 2 -x* 2 



z:=a -l/(l-B) 



W(Z 



^ a 4-n 2 x* 2 ' 

(C"|x), (4.1) 

\ 2 — n A — n z J 



Then equation ( |2.8|) takes the form 

w(z) 2 + w(-z) 2 + nw(z)w(-z) = [i sin 2 n6 + a 2e/{1 ~ e) z 2 /(2 - n) . (4.2) 



In the continuum limit a — » 0, the second term on the r.h.s. of ( |4.2| ) can be neglected 
provided h < 9 < 1, or0<n<2. Then the solution of ( [4.2|) can is written in parametric 
form as 

z = M cosh r , 
w(z) = — M l ~ cosh(l — 9)t , <U> 

where M = C / u 1 / 2< - 1_6 '- ) . The value of the constant C is fixed by the normalization of the 
disk partition function $. One possible choice is 

M e 

= — cosher, a 2 $ = -w(z), M 2 ~ 2e = (4.4) 
9 



5 Note that in the dense phase the boundary has anomalous dimension. This is a consequence 
of the fractal structure of the boundary in this phase. 
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As a function of z, the loop amplitude w(z) has a branch cut along the interval [—00, — M\. 
The solution (|4.3|) corresponds to Liouville gravity with matter central charge 



1-6 



1-6 



(4.5) 



The susceptibility which is by definition the partition function of the loop gas on a 

dynamically triangulated sphere with two punctures, scales as u ~ n~ lst \ with 



6 
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1-6 



(4.6) 



4-2. Boundary two-point functions of the L-leg operators 



The continuum limit of the boundary two-point functions is obtained as a triple scaling 
limit in x, ( and £, in which the area of the triangulation as well as the lengths of the 
Neumann and JS boundaries diverge. The point x = x*, ( = £*, C = C* is a singular point 
of equation (|3.20| ), where the r.h.s. vanishes: 



w(o + -w(o + -(w(n - a = . 

y y 



(4.7) 



We define the renormalized cosmological constant z and loop amplitude with JS boundary 
conditions as 



z : = a 



•1/(1-0) 



c-c 



w(~z):=a- 1 (W(C)-W(C)) 



In the continuum limit a — > 0, the linear terms in z and z can be neglected, since they are 
multiplied by a 6 ^ 1-0 ^, and functional equation Q3.20|) takes the form 



DHz)D^(-z) 



-w(z) w(—z) w(z). 

y y 



(4.9) 



Together with ( |3.18| ), this equation yields a linear relation between and Dq. 



D^(z) 



w(z) H — w(— z) H — w(z) 

y y 



m-z) 



(4.10) 



Equation ( |4.10| ) is the key result of this paper. It allows to determine, up to a normalization, 
the boundary two-point functions with D± and Dq. The two-point functions with L > 1 
are then easily evaluated from the recurrence equations Q3.8Q , which have, in the continuum 
limit, the form 



D L (z + iQ)-D L (z-iO) = [w(z + iO)-w(z-iQ)]D L - 1 (z) (z < -M) . (4.11) 



14 



The solution, as a function of fi, z and z, is expected to be of the form 

D L (jji,z,z)=n 1 -^{y/ii) 2 ^- 2 D L (z/M,z/M), (4.12) 

where the exponent 7 s t r is given by Q4.6|) , A£ is the boundary gravitational dimension of 
the L-leg operator, and M = (2^) 1 / (2_2e) . 

Before giving the complete solution, we are going to solve a simpler problem: to 
determine the scaling behavior of the two-point function when \x = z = 0. Restricted in 
this way, the two-point function ( f4.12|) reduces to a power of the only non-zero coupling z: 

Di ^(l-0)(2Af- 7str ) (413) 

From here one can determine the conformal dimensions of the L-leg operators using the 
KPZ map Q. 



4-3. Evaluation of the L-leg critical exponents for N-JS boundary conditions 

At fx = z = 0, the solution for the observables Dl and w must be given by powers of z, 

Di~z aL , D ] l^z pL , w~z l -°. (4.14) 

Then ( fl.9p yields three identities relating (3q, [3\ and y. The first one follows from comparing 
the powers of z on both sides: 

P 1 + ao = l-9. (4.15) 

The other two identities arise when equating the imaginary part of both sides of (|4.9| ) 
respectively for z > and z < 0: 

sin(7rao) = sin#7r, sin(7r/3i) = — sin(^7r). (4-16) 

y 

These equations determine y as a function of a®: 

shi7r/3i sin7r(a + 0) . ^ 

y = = : • (4.17) 

sin noto sin nao 



We see that the expression Q4.17|) for y is of the form (|1.6| ) with 

a = r6. (4.18) 

To determine ccq, we must invert the multi-value function y(ao) = y(«o + !)• The 
relevant branch is given by the lowest positive value of «o, which is must be in the interval 

0<ct o <l or 0<r<^. (4.19) 

V 
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The fact that ao is positive is a consequence of the representation (|3.19|) of Dq. The 
singularity of the l.h.s. cannot not be weaker than the singularity of each of the terms, 
hence (3q < ao. Since ao + (3q = 0, this means that «o is positive and (3q is negative.! 

Once we determined the critical exponent for L = 0, the other exponents follow from 
the recurrence relation ( f4.11| ): 

a L = L(l-9)+r9 j3 L = L(l - 9) - r9 . (4.20) 

Comparing these values with ( f4.13| ) we find for the gravitational dimensions and A^; 

± _ a L -9 _ L(l -9) + r6-9 _ 

(1.21) 

A —j. —<f — £, 



A 
A 



L 2(1 - 9) 2(1 - 9) 

(3 L -9 L(l-9)-r9-9 



2(1-9) 



2(1-9) 



By the KPZ scaling relation ( |1.1[ ), and using the symmetry /i_ r _ s = h rs , we find the 
conformal weights of the L-leg operators with mixed Neumann-JS boundary conditions: 



(4.22) 



in accord with 18 



When r = 1, then y = n and the loops touching the boundary have the same fugacity 
as the loops in the bulk. In this case one obtains the well known L-leg exponents with 
Neumann- Neumann boundary condition on flat [19| and dynamical |22] lattices. 



h "L* = Vl/2+L 



4-4- Complete solution and relation to Boundary Liouville theory 



In the previous subsection we determined the prefactor in (|4.12| ). Now we will evaluate 
the scaling functions Dj_. First, using the expression (|4.3f ) of the loop amplitude w(z) and 
the identity 



[w(t) + y~ l w(T + Z7r)]/M 1_e = -C cosh[(l - 9)t + inr9] 
with C = smir9/ sin(r + 1)tt9, we write ( [4.10| ) as 



(4.23) 



Dll(r) 



C cosh[(l-^)r + Z7rr^] + 



to u 



yM 1 - 



D]}(T + m). 



(4.24) 



6 Let us stress that this argument is justified only for non-negative couplings, when all terms 
in the series are non-negative. 
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Next, we change the variable z — > r so that 

cosh(l-0)f 



After that equation ( |4.11|) takes the form 



sin tt9 
sin nrO 



cosh(l — 0)t 



t)\\ ( r , f) = -C [cosh ((1 - 6)t + mr9) + cosh(l - 0)?)] DH (r + iir, r) . 



(4.25) 



(4.26) 



Let us remark that equation (|4.25|) is just a change of variable and does not involve any 
assumption about the the boundary one-point function with JS boundary conditions w(z), 
because the function z(f) is not yet determined. On the other hand, the only solution 
compatible with the world-sheet CFT is z ~ Mcoshf. Indeed, there is only one non-zero 
boundary one-point function, that of the Liouville-dressed identity operator, which is given 

by©- 

For generic r, one can recognize in (|4.26| ) the functional equation derived by V. Fa- 
teev, A. Zamolodchikov and Al. Zamolodchikov for the boundary two-point function in 
boundary Liouville theory jl5] . This equation also appeared as the first member of an infi- 
nite series of functional identities characterizing the boundary ground ring in 2D quantum 
gravity JZ§[fZp]. 

The boundary two-point function in Liouville gravity depends on the the target-space 
momentum p and the two boundary parameters r and f .0 It is given, up to a normalization 
factor that depends only on p, by |15[ 



D(p;r,f)=^ 2 D(p;r,f), 



D(p; r,r) = exp 



dt 

T 



sinh(7rpt/6 2 ) cos(rt) cos(rt) p 
sinh(7rt) sinh(7rt/6 2 ) nt 



(4.27) 



The scaling function D(p; r, f) satisfies the identity 

D(p + b 2 ; r, f) = \ [cosh (6 2 r t inp) + cosh (6 2 r)] D(p; r ± ztt, r) , (4.28) 
which is the same as ( [4.26] ), with b 2 = 1 — 9, p = —r6 and C = 1/2. Therefore 

D^(T,f) = D(-er;T,f), D l }(T,T)=b(-er + l-e-,T,T) (4.29) 
is a solution of (|4.26|) . The relation (|3.18j ) then implies 



(4.30) 



7 Our notations are related to the notations of pl| by 1/6 + 6—2/3 = p/b, r = ns/b, b 2 = 1 — 9. 
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Is this the physical solution that corresponds to the series expansion ( |3.4p ? Equation 



Q4.28| ) does not determine uniquely the two-point function. In order to specify the unique 



solution, FZZ |TH] used the duality symmetry, which supplies another equation of the 



same type, as well as the symmetry in r <-> f. Neither of these symmetries is satisfied in 
the microscopic theory. If we assume that the physical solution enjoys these symmetries 



in the scaling limit, then (|4.29|) is the unique solution. In the particular cases r = 1 and 
r = (l — 9)/29, the function ( |4.3U| ) reproduces correctly the expressions obtained previously 



for the disk partition function with mixed Neumann/ Neumann and Neumann/Dirichlet 
boundary conditions. 

Once the two-point functions with L = are known, the rest can be determined from 
the recurrence equation ( |4.11| ), which can be cast into the form 

b L {r + m,f) - D L (r -in,?) = -2zsin7r# sinh[(l - 6)t] D L -i(t, f). (4.31) 

This equation is compatible (up to normalization) with G4.28|) upon the identification 

p\ = -dr + (1 - 9)L and pi = dr + (1 - 9)L , (4.32) 

respectively for and Dj^. When L > 1, the operators and E>i have different 
conformal weights. The two-point functions Dq and Dq are related by Liouville reflection, 
and correspond to the 'physical' and 'unphysical' Liouville dress of the same boundary 
matter field with conformal weight h r ^ r . 

4-5. Dirichlet boundary conditions and twist operators 

The Dirichlet boundary conditions for the 0(n) model is defined by fixing the 0(n) spin 
on the boundary to point to given direction, say S = (1, 0, ... 0). This is a particular case 
of the JS boundary condition, obtained by taking y = 1, or equivalently 

r=±^. (4.33) 

This case deserves special attention, because here we can compare the general formula 
(ODD with the exact results found in |E3 IS3] Eg . 




Fig. 7: A loop gas configuration for Dirichlet boundary condition. 
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The Dirichlet boundary condition for the 0(n) model leads to the same loop gas 
expansion as the Neumann boundary condition for the SOS model, which was first studied 
22] and then given a world sheet CFT interpretation in p4| and |23|. Namely, it is 



m 



assumed that each point of the boundary is an endpoint of an open line, as sketched in 
Fig. 7. The boundary condition changing operator T separating Dirichlet and Neumann 



boundary conditions was called in |24j] twist operator, by analogy with the gaussian field. 
The correlation function of two twist operators, 



Cl(z, z) 



N 



disk 



(4.34) 



was first evaluated in [02 1, eqs.(4.34)-(4.37) of that paper. Afterwards this solution was 
identified [f21| as special case of the FZZ two-point function ( 4.27 ) with p = (1 — #)/2, 
which corresponds to r given by ( |4.33|) . 

It is easy to see that the sum over the configurations with open lines as the one in Fig. 
7 can be interpreted, for this particular value of r, either as the loop expansion for Dq, 
or as the loop expansion for d\ . Indeed, if we connect pairwise the endpoints of the open 
lines, as shown in Fig. 8a, we obtain a configuration of the loop expansion for Dq . The 
Boltzmann weights also match under the condition that all loops that touch the boundary 
have fugacity y = 1. Alternatively, we can leave the first and the last open line endpoint 
free and connect the rest of the endpoints pairwise, as is shown in Fig. 8b. Then the first 
and the last points are connected by an open line, and we obtain a configuration of the 
loop expansion for Therefore, even at microscopic level, 



= = (y = l). 





a b 
Fig. 8: The two ways of closing the open lines. 



The two-point function of the boundary twist operator, 0(r, r), is obtained as the 
solution of a quadratic functional identity, eq. (4.25) of |P2"| , which is identical to ( fO| ) with 
y = 1. In this particular case the solution can be found without additional assumptions 
about symmetry, and the result |22| coincides with (|4.27|) for p = (1 — 0)/2. 



19 



The conformal weights of the excited twist operators Tc, or the boundary L-leg op- 
erators with mixed Neumann/Dirichlet boundary conditions, were identified in |23| as§ 

hr L = h 0jL+1/2 ■ (4.35) 

On can check that Ht l = h^l +1 = hj^ , with r given by Q133|) . The different Kac-table like 
identifications of these operators are possible because of the ambiguities of the represen- 
tation flTg) . 

Thus the results obtained here for general r are in full agreement with those of 
22,^4,23]. There is however a difference in the interpretation of the results, which is 



due to the different form of the microscopic loop equations. Compared to equation (4.18) 
of 1 2^1 , our equation ( |3.20| ) contains an extra term W, which takes into account the self- 



touchings of the JS boundary. This term was omitted in ||22j| . hence the self-touchings of 
the Dirichlet boundary were not taken into account there. As a result, the two-point func- 
tion ( |4.34| ) , evaluated in |22 ] , describes a sum over surfaces with self-touchings allowed for 



the Dirichlet boundary and forbidden for the Neumann boundary. This explains the puz- 
zling observation, made in ||22j| , that the fractal dimensions of the Dirichlet and Neumann 
boundaries are different for 9 ^ 0. 

We find here more natural to define the sum over surfaces so that both segments of the 
boundary have the same dimension, which is the case when the contact term in question 
is taken into account. 



5. Conclusions 

In this paper we evaluated the boundary two-point function for the 0{n) loop model on 
the dynamically triangulated disc with presumably the most general boundary conditions, 



constructed recently in |18|. We restricted ourselves to the dense phase of the loop gas, 
where both the bulk and the boundary are critical and the only parameters of the theory 
are the bulk and boundary cosmological constants. The scaling behavior of the two-point 
function confirms the L-leg exponents ( f4.22| ) conjectured in |[18|| . Our result for the two- 
point function implies the symmetry 

Si <- Si, y <- n-y. (5.1) 

which looks quite natural given the definition ( |3.2[ ) of these operators and resembles the 
duality symmetry that exchanges Dirichlet and Neumann boundary conditions. In the 



In [23 1, the conformal weights for the excited twist operators were actually evaluated for the 



dilute phase, where they are given by h L+1 / 2 ,o- 
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parametrization (|TTE|), exchanging y and n — y is equivalent to changing the sign of r, due 
to the identity 

y{r) + y(—r) = n. (5-2) 
In the loop gas formulation, the symmetry (|5.1|) is spelled out as 

{blobbed, r} «-> {unblobbed, — r} (5-3) 

and is respected by the exponents (|1.5|) . 



The dilute phase is more intricate because for each y there is a one-parameter family 
of boundary conditions and a fine tuning of the matter coupling constants should be done 
bth in the bulk and on the boundary. Our preliminary results for the scaling dimensions 
30(1 seem to be compatible with the unpublished results of Jacobsen and Saleur. 

Finally, let us mention that the 0(n) model coupled to 2D gravity can be viewed as 
solvable model of bosonic string theory with curved target space, representing the (n — 1)- 
dimensional sphere. The continuous spectrum of D-branes in this theory is presumably 
related to the fact that the target space curvature, 

R = (n-l)(ra-2), 

is negative in the interval 1 < n < 2. We believe that all our results can be also obtained 
also on the basis of the dual 0{n) invariant matrix model ||. 

Acknowledgements 

The author thanks J. Jacobsen and H. Saleur for discussions and explaining me some of their 
unpublished work on this subject. Part of this work was carried out at the Kavli Institute for 
Theoretical Physics, Santa Barbara, supported in part by the National Science Foundation under 
Grant No. PHY99-07949, in the course of the program on Stochastic Geometry and Field Theory: 
From Growth Phenomena to Disordered Systems. This work has been partially supported by the 
European Network ENRAGE (contract MRTN-CT-2004-005616) and the ANR program GIMP 
(contract ANR-05-BLAN-0029-01). 



21 



References 



[1] V. K. Kazakov, "Percolation on a fractal with the statistics of planar Feynman graphs: 

Exact solution", Mod. Phys. Lett. A 4, 1691 (1989). 
[2] V. A. Kazakov, "Exact solution of the Ising model on a random two-dimensional 

lattice", JETP Lett. 44, 133 (1986) [Pisma Zh. Eksp. Teor. Fiz. 44, 105 (1986) ]. 
[3] D. V. Boulatov and V. A. Kazakov, "The Ising Model On Random Planar Lattice: 

The Structure Of Phase Transition And The Exact Critical Exponents," Phys. Lett. 

186B, 379 (1987). 

[4] B. Duplantier and I. Kostov, Phys. Rev. Lett 61 (1988), 1433; Nucl. Phys. B 340 
(1990) 491. 

[5] I. K. Kostov, "O(n) Vector Model On A Planar Random Lattice: Spectrum Of Anoma- 
lous Dimensions," Mod. Phys. Lett. A 4, 217 (1989). 

[6] I. K. Kostov, "The ADE face models on a fluctuating planar lattice" Nucl. Phys. B 
326, 583 (1989). 

[7] I. K. Kostov, "Strings with discrete target space," Nucl. Phys. B 376, 539 (1992) 
|arXiv:hep-th/9112059| . 



[8] A. Polyakov, "Quantum Geometry Of Bosonic Strings" , Physics Letters B103 (1981), 
207. 

[9] V. Knizhnik, A. Polyakov and A. Zamolodchikov, Mod. Phys. Lett. A3, 819 (1988); 
F. David, Mod. Phys. Lett. A3, 1651 (1988); J. Distler and H. Kawai, Nucl. Phys. 
B321, 509 (1989). 

[10] F. David, "Conformal field theories coupled to 2-D gravity in the conformal gauge", 
Mod. Phys. Lett. A 3, 1651 (1988); J. Distler and H. Kawai, "Conformal field theory 
and 2-D quantum gravity or who's afraid of Joseph Liouville?", Nucl. Phys. B 321, 
509 (1989). 

[11] J. L. Cardy, "Conformal Invariance And Surface Critical Behavior," Nucl. Phys. B 
240, 514 (1984). 

[12] B. Duplantier, "Conformal fractal geometry and boundary quantum gravity," ||arXiv:math- 
| ph/030303l . 

[13] B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982); J. Stat. Phys. 34, 731 (1984). 

[14] See, for example, G.F. Lawler, Conformally Invariant Processes in the Plane, Mathe- 
matical Surveys and Monographs, Volume 114 (American Math. Soc, April 2005). 

[15] V. Fateev, A. B. Zamolodchikov and A. B. Zamolodchikov, "Boundary Liouville field 
theory. I: Boundary state and boundary two-point function", ||arXiv:hep-th/0001012 |. 

[16] B. Ponsot, J. Teschner, "Boundary Liouville field theory: Boundary three point func- 
tion" , Nucl. Phys. B 622 (2002) 309, ||arXiv:hep-th/0110244j . 



22 



[17] 

[18] 
[19] 
[20] 
[21] 
[22] 

[23] 

[24] 

[25] 



[26] 
[27] 

: 

[28] 
[29] 
[30] 



K. Hosomichi, "Bulk-Boundary Propagator in Liouville Theory on a Disc", JHEP 
0111, 044 (2001), ||arXiv:hep-th/0108093| . 

J. Jacobsen and H. Saleur, "Conformal boundary loop models" , [|arXiv:math-ph/0611078 
B. Duplantier and H. Saleur, Phys. Rev. Lett. 58, 2325 (1987). 

H. Saleur and M. Bauer, Nucl. Phys. B 320, 591 (1989). 

J. L. Cardy, "The O(n) model on the annulus," ||arXiv: math-ph / 0604043) . 

V. Kazakov and I. Kostov , "Loop gas model for Open Strings", Nucl. Phys. B 386, 
520 (1992). 

I. K. Kostov, B. Ponsot and D. Serban, "Boundary Liouville theory and 2D quantum 
gravity", Nucl. Phys. B 683, 309 (2004), [|arXiv:hep-th/0307189|| . 

I. K. Kostov, "Boundary correlators in 2D quantum gravity: Liouville versus discrete 
approach," Nucl. Phys. B 658, 397 (2003) ||arXiv:hep-th/0212194j . 

A. Nichols, V. Rittenberg and J. de Gier, J. Stat. Mech. 0503, P003 (2005), 
arXiv : cond-mat / 04 1 1 5 1 2| ; A. Nichols, J. Stat. Mech. 0601, P003 (2006), Ihepl 



th/0509069| ; A. Nichols, J. Stat. Mech. 0602, L004 (2006), [|arXiv:hep-th/0512273|| . 
P. A. Pearce, J. Rasmussen and J. B. Zuber, "Logarithmic minimal models," J. Stat. 
Mech. 0611, P017 (2006) ||arXiv:hep-th/0607232 l. 

A. Y. Alekseev and V. Schomerus, Phys. Rev. D 60, 061901 (1999); ||arXiv:hep 
th/9812193|| . 

M. Bershadsky and D. Kutasov, "Scattering of open and closed strings in (1+1)- 
dimensions", Nucl. Phys. B 382,213, ( 1992), |hep-th/ 9204049 . 

I. K. Kostov, "Boundary ground ring in 2D string theory", Nucl. Phys. B 689, 3 
(2004), [|arXiv:hep-th/0312301|] . 

I. Kostov and Al. Zamolodchikov, work in progress. 



23 



